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ABSTRACT 

We present wide-field photometry of resolved stars in the nearby Sculptor dwarf spheroidal galaxy using 
CTIO/MOSAIC, going down to the oldest Main Sequence Turn-Off. The accurately flux calibrated wide 
fieldColour-Magnitude Diagrams can be used to constrain the ages of different stellar populations, and also 
their spatial distribution. The Sculptor dSph contains a predominantly ancient stellar population (>10 Gyr 
old) which can be easily resolved into individual stars. A galaxy dominated by an old population provides 
a clear view of ancient processes of galaxy formation unimpeded by overlying younger populations. By 
using spectroscopic metallicities of RGB stars in combination with our deep Main Sequence Turn-Off 
photometry we can constrain the ages of different stellar populations with particular accuracy. We find that 
the known metallicity gradient in Sculptor is well matched to an age gradient. This is the first time that 
this link with age has been directly quantified. This gradient has been previously observed as a variation 
in Horizontal Branch properties and is now confirmed to exist for Main Sequence Turn-Offs as well. It is 
likely the Sculptor dSph first formed an extended metal-poor population at the oldest times, and subsequent 
more metal-rich, younger stars were formed more towards the centre until the gas was depleted or lost 
roughly 7 Gyr ago. The fact that these clear radial gradients have been preserved up to the present day is 
consistent with the apparent lack of signs of recent tidal interactions. 

Key words. Galaxies: dwarf - Galaxies: evolution - Galaxies: stellar content - Galaxies: Local Group - 
Stars: C-M diagrams 
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1. Introduction 

The Sculptor dwarf spheroidal galaxy is a 
faint (My^-11.2), well studied system in the 
Local Grou p. It has a tidal radius of 76.5 arcmin 
on the sky dlrwin & Hat zidimitrioul [l995h . which 
corre sponds to 1.9 kpc at a distance of 86+5 
kpc dPietrzvnski et all UW^i . It is located at high 
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galactic latitude (b=-83°) wit h a systemic ve locit 
of Whei = +1 10.6 ± 0.5 km/s dBattaglia et alJl200 
and suffers from rela tively low amounts o f redden- 
ing, E(B-V)=0.018 dSchlegel et al.lll998h . Sculptor 
was the first early-typ e dwarf galaxy discovered 
around the Milky Way (Shap leylll938l) . Since then 
it has been the target of numerous studies of its 
resolved stellar populations. 

Early work on the radial a nd two-dimensional 
structure of the S cul ptor dSph (|Demers et al.l|l980t 
lEskridgd 119888111 llrwin & Hatzidimitrioul 1 19951) 
uncovered a complex spatial structure, with a 
radially increasing eccentricity. Furthermore, the 
presence of a spatial asymmetry was suggested in 
the East/West direction in observed star counts, 
of the order 30% above the level expected for a 
symmetric profile (Eskridge 1988b ). Using Colour- 
Magni tude Diagrams (CMP). iNorris & Besselll 
(1 19781). interpreting d a ta fro m iHodgel (1 19651) 
and iKunkel & Demersl dl977l) . argued that the 
Red Giant Branch (RGB) is wider than can be 
explained by photometric uncertainties, suggesting 
the presence of internal age and/or abundance vari- 
ations. The first CMDs using CCDs in a field just 
outside the core radius confirmed the large RGB 
spread and extended down to the Main Sequence 
Turn-Off (MST O) region, deter mining an age range 
of 13+2 Gyr (iDa Costal 11984). A very deep HST 
CMD of a small field of view (~2') well outside 
the centre of Sculptor, going down 3 magnitudes 
below the oldest MSTO, accurately confirmed the 
ancient age (15+2 Gyr) of the bulk of the stars in 
Sculp tor. This is similar to th e age of globular clus- 
ters dMonkiewicz et al.lll999l) . From the same HST 
data set it was shown using CMD synthesis analysis 
that Sculpt or disp l ays an extended Star Formation 
History (Dolphin 2002), with most of the stars 
having formed 8-15 Gyr ago, but with a small (and 
highly uncertain) tail reaching down to more recent 
times. Ground based wide-field imaging studies, 
covering a larger fraction (a;30' x 30') of Sculptor, 
including the centre, found that the Horizontal 
Branch (HB) morphology chang es significantly 
with distance from the cent re dMaiewski et al.l 
1999: lHurlev-Kelleret aljfl999h . 
On the RGB the presence of the age-metallicity 
degeneracy is a problem for the accurate analy- 
sis of the properties of the Sculptor dSph. The 
age-metallicity degeneracy on the RGB results in 
two populations displaying nearly identical RGB, 
if the second population has an age a factor of 3 
higher and a metallicity a factor of 2 lower than 
that of population 1. Therefore, spectroscopic 
observations are needed to remove the degeneracy 



on the RGB and allow an accurate analysis of 
Sculptor properties. 

Wide-field medium resolution Can triplet spec- 
troscopy of RGB stars in Sculptor have indepen- 
dently confirme d the presence of two distinct stellar 



comp o nents dTolstoy et al.l I2004t iColeman et alj 
120051: IWestfall et alJ 120061) . with different spa- 
tial distributions, kinematics and metallicities. 
There is a spatially extended component which 
is more metal-poor (-2.8<[Fe/H]< -1.7), and 
a more concentrated, metal-rich component 
(-1.7<[Fe/H]< -0.9) which also has a lower 
velocity dispers i on than the metal- poor component 
(Battaglia 2007; Battaglia et al. 2008). 
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Fig. 1: Coverage of the photometric data across the 
Sculptor dwarf spheroidal. Black squares (no. 1-5,9-12) 
denote the CTIO 4m fields observed. The 15 smaller 
blue squares show the CTIO 0.9m calibration fields that 
were observed. The green dashe d ellipse is the tidal ra- 
dius o f Sculptor, as determined by Irwin & Hatzid imitrioul 



High resolution spectroscopic studies in Sculptor 
have been carried out to determine detailed abun- 
dances of a range of chemical ele ments for small 
numb ers of indi vidual RGB stars dShetrone et all 
120031 [5 stars]; iGeisler et al l 120051 [4 stars]) us- 
ing VLT/UVES and more recently with wide-field 
multi-fibre spectrograph VLT/FLAMES in high res- 
olution mode (Hill et al., in prep; see iTolstov et al.l 
(2009)) for a sample of 93 stars. These studies 
are able to constrain the chemical evolution of 
Sculptor. This can then be compared to the other 
Local Group galaxies, including the Milky Way. 
For example, the [a/Fe] abundances are sensitive 
to the relative rates of SNe II and SNe la and can 
thus be used to trace the star formation timescale. 
Thus, the [Fe/H] and [ff/Fe] abundances obtained 
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from spectroscopic studies can be combined with 
deep photometry to accurately determine the ages 
of the Sculptor stellar populations, and the relation 
between Star Formation Rate and the chemical en- 
richment processes. Obtaining accurate time scales 
for star formation and the build up of abundance 
patterns in Sculptor also allows us to determine the 
early history of Sculptor, and also until which time 
Sculptor could contribute to the build-up of the 
Milky Way. 

In this work we present deep wide-field carefully 
calibrated photometry of the Sculptor dSph galaxy, 
earned out with MOSAIC on the CTIO 4m/Blanco 
telescope. These CMDs reach down to the oldest 
MSTO, for an area covering «80% of the tidal ra- 
dius of Sculptor. Accurately calibrated photometry 
is essential in determining ages, since the position 
of the MSTO directly correlates with the age of the 
stellar population. From deep MSTO photometry 
over a large region we can for the first time study 
the effect of age on the physical properties of the 
galaxy. Combining the results of spectroscopic 
surveys with our MSTO photometry also allows us 
to measure the age-metallicity relation in Sculptor 
dSph galaxy. 

The paper is structured as follows: in section [2] 
we present the observations and data reduction. In 
section [3] we describe how we obtained and cali- 
brated the photometry and the resulting structural 
properties. Section |4] describes the analysis and in- 
terpretation of different evolutionary features in the 
CMDs. Finally, section [5] discusses the conclusions 
and their implications in terms of galaxy formation. 
In a subsequent paper we will present the detailed 
Star Formation History analysis of the Sculptor 
dSph. 

2. Observations & Data Reduction 

Deep optical photometry of the Sculptor dSph in the 
B, V and I bands was obtained using the CTIO 4- 
m MOSAIC II camera over 10 nights in September 
2008 and November 2009. In order to ensure accu- 
rate photometric calibration without using too much 
4m telescope time we obtained service mode obser- 
vations with the 0.9m CTIO telescope over 3 photo- 
metric nights. O bservations we r e mad e of Landolt 
standard fields (iLandolll 120071 [l99l covering a 
range of different airmass and colour. 
The positions of fields observed with the 4m (big 
black squares) and 0.9m (small blue squares) tele- 
scopes are shown in Fig. [T] relative to the Sculptor 
centre. B,V and I photometry was obtained for the 
central 5 pointings (1,2,3,4,5), whereas the outer 



pointings (9,10,1 1,12) were only observed in V and 
I bands. The observation logs for Sculptor obser- 
vations are given in Tables Q] and [2] The observa- 
tion logs for standard star calibration fields with 
both the 4m and the 0.9m telescopes are given in 
AppendixlAl 

2.1. CTIO 4m MOSAIC 

The CTIO MOSAIC II camera has an array of eight 
2048x4096 pixel CCDs, each of which can be read 
out through one or two amplifiers. The CCDs com- 
bine to form a 8192x8192 pixel image correspond- 
ing to a field of view of 36' x 36' on the sky. For 
our observing runs, one amplifier on one CCD was 
found to be inoperable, such that we were forced to 
read out MOSAIC II with 8 amplifiers only. This 
resulted in a readout overhead of 2m40s per im- 
age. An observing strategy was chosen involving 
non-dithered long (600s) exposures to be stacked to- 
gether. We chose not to dither images, in order to be 
able to directly co-add images to obtain the deepest 
photometry possible. Furthermore, not dithering al- 
lowed us to apply accurate position dependent aper- 
ture corrections for each pointing. The loss of area 
incurred by not dithering exposures is roughly 2 per- 
cent, which means it will not materially affect the 
conclusions in this paper. However, the very inner 
part of Sculptor (r e ;; <0.05 deg) does suffers from 
incomplete coverage due to the presence of CCD 
gaps there. Since the average saturation level of the 
CCDs is ^43000 counts (corresponding to Ba;18, 
V»T7.5 and 1*17.5 for 600s integration) additional 
shorter exposures (90s + 10s) were taken to sample 
the bright stars in each field. 

For the data reduction we used the data 
reduction program, wh ich includes the MSCRED 
package ( Valdesl l 19981) . designed for working with 
mosaics of CCD's. The standard data reduction 
steps include correcting for bias and Hatfield, re- 
moving bad pixels and obtaining precise astrome- 
try. For wide-field images there are a number of 
ext ra steps whic h needed to be taken, as outlined 
in (Valdes 2002). 

The 8 chips of each MOSAIC image suffer from 
electronic crosstalk which causes ghost images. 
Each pair of CCD's are controlled by one CTIO 
Arcon controller, which causes one CCD to show 



1 IRAF is distributed by the National Optical 
Astronomy Observatories,which are operated by the 
Association of Universities for Research in Astronomy, 
Inc., under cooperative agreement with the National 
Science Foundation. 
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Table 1 : List of fields observed in the Scl dSph with the 4m CTIO Blanco telescope. 



Date 


rield 


KA 




Filter 


exp time 


seeing 


airmass 






J2000 


J2000 




sec 






"\ i \ / \l \ XT 1 t \ 

2009 Nov 19 


Sell 


01:00:03.96 


-33:41:30.48 


B 


1800, 90, 10 


1.0-1.3 


1.00-1.01 


2009 Nov 19 








V 


1800, 90, 10 


1.1-1.2 


1.01-1.02 


2009 Nov 19 








I 


2400, 90, 10 


1.0-1.3 


1.03-1.07 


2008 Sep 05 


Scl2 


r» 1 ci r\o 

01:02:52.08 


1 At i~\ A r\c\ 

-33:41:24.00 


B 


1500, 60 


1.6-2.2 


1.00-1.06 


2009 Nov 20 








B 


2400 


1.6-2.2 


1.00-1.06 


2008 Sep 05 








V 


1500, 60, 30 


1.0-1.1 


1.00-1.01 


2009 Nov 21 








V 


1800 


1.0-1.1 


1.00-1.01 


2008 Sep 05 








I 


a f\f\ f f\ 1 f\ 

2400, 60,10 


1.2-1.3 


1 AA 1 C\^l 

1.00-1.07 


2009 Nov 23 


Scl3 


00:57:15.48 


-33:41:22.92 


B 


900, 60 


1.2-1.4 


1.62-1.70 


2009 Nov 23 








B 


600 


1.2-1.4 


1.62-1.70 


2008 Sep 06 








V 


900, 60, 30 


1.0-1.2 


1.01-1.02 


-\ / \ / W \ XT 1 

2009 Nov 21 








V 


1800 


1.0-1.2 


1.01-1.02 


2008 Sep 06 








I 


1800, 60, 10 


0.9-1.0 


1.03-1.08 


2009 Nov 21 








I 


2400 


0.9-1.0 


1.03-1.08 


2008 Sep 06 


Scl4 


01:00:03.96 


-33:06:30.60 


B 


1 c r\r\ r r\ 1 a 

1500, 60, 10 


0.9-1.0 


1.01-1.03 


2008 Sep 06 








V 


900, 60, 10 


1.0-1.5 


1.10-1.30 


2009 Nov 21 








V 


1800 


1.0-1.5 


1.10-1.30 










T 
1 




1 n 1 ^ 




2008 Sep 06 


Scl5 


01:00:03.96 


-34:16:30.36 


B 


900, 60, 30 


0.8-1.1 


1.00-1.01 


2008 Sep 06 








V 


900,60, 30 


0.7-1.1 


1.00-1.01 


2009 Nov 22 








V 


1800 


0.7-1.1 


1.00-1.01 


2009 Nov 22 








I 


2400, 90, 10 


1.1-1.3 


1.01-1.04 


2009 Nov 22 


Scl9 


00:57:16.92 


-33:06:23.40 


V 


1800, 90, 10 


1.1-1.2 


1.06-1.11 


2009 Nov 22 








I 


2400, 90, 10 


1.0-1.2 


1.14-1.28 


2009 Nov 23 


ScllO 


00:57:14.40 


-34:16:22.80 


V 


1800, 90, 10 


1.1-1.2 


1.10-1.17 


2009 Nov 23 








I 


2400, 90, 10 


1.0-1.2 


1.21-1.39 


2009 Nov 23 


Sclll 


01:02:53.16 


-34:16:23.52 


V 


1800, 90, 10 


0.9-1.0 


1.00-1.02 


2009 Nov 23 








I 


2400, 90, 10 


0.8-0.9 


1.00-1.02 


2009 Nov 22 


Scl 12 


01:02:51.00 


-33:06:24.12 


V 


1800, 90, 10 


1.0-1.1 


1.29-1.43 


2009 Nov 23 








I 


2400, 90, 10 


0.9-1.0 


1.03-1.07 



ghost images from stars in a different CCD that are 
saturated or nearly saturated. This crosstalk was re- 
moved using correction terms provided by NOAO. 
The flat-field obtained for the MOSAIC camera is 
made up of different chips. Instead of normaliz- 
ing each chip individually, all chips were normal- 
ized using the average mean across the whole im- 
age. The normalization factor for each chip then de- 
pends on its exact gain, and can be used to correct 
the science images both for the variation of sensi- 
tivity across each chip and to normalise the gain be- 
tween different chips to a common value, thus ob- 
taining a final image which has one uniform aver- 
age gain. The standard flat-field obtained using the 
dome lamp was improved with the use of sky flat- 
fields. Sky flat-fields better approximate uniform il- 
lumination than possible with the dome lamp, as 
light enters the telescope in the same way as dur- 
ing science observations. The sky flat-field provides 



a correction on the dome flat-field, creating a mas- 
ter flat-field which is used to correct the science im- 
ages. The MOSAIC I-band images also suffer from 
fringing, caused by constructive and destructive in- 
terference of night sky lines reflected in the CCD 
substrate. These fringe patterns were removed us- 
ing an archival I-band fringe frame. 
We created a bad pixel mask to flag saturated and 
negative pixels, cosmic rays and any damaged ar- 
eas on each CCD. Cosmic rays were detected with a 
moving 7x7 pixel block average filter with the cen- 
tral target pixel excluded. The average value of the 
block is a prediction for the value of the target pixel, 
and the median of a background annulus around the 
block is used to protect cores of stars from being 
flagged as bad pixels. If the difference between the 
average block value and the central pixel value was 
larger than a threshold, the target pixel was flagged 
as a cosmic ray and added to the bad pixel mask. 
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Table 2: List of fields observed in the Scl dSph with the 0.9m CTIO telescope. 



Date 


Field 


RA 


DEC 


Filter 


exp time 


seeing 


airmass 






TO AAA 
J2000 


TO AAA 
J2000 




sec 


n 




2008 Oct 3 1 


Sell N 


n 1 f\r\ no n s~ 

01:00:03.96 


00 on a a 

-33:32:30.48 


B 


300 


1.5 


1.02-1.05 


2008 Oct 3 1 








V 


300 


1.4 


1.02-1.04 


2008 Oct 3 1 








I 


600 


1.3 


1 nn 1 no 

1.00-1.03 


r\l~\C) XT AO 

2008 Nov 02 


Sell S 


a 1 nn no n/" 

01:00:03.96 


00 c n on ao 

-33:50:30.48 


B 


300 


1.5 


1.06-1.09 


2008 Nov 02 








V 


300 


1.3 


1.05-1.08 


nAAO XT AO 

2008 Nov 02 








I 


600 


1.3 


1.04-1.07 


2008 Oct 3 1 


Scl2 N 


ni no c — i nn 

01:02:52.08 


00 a nn 

-33:32:24.00 


B 


300 


1.7 


1.01-1.03 


in/Mi /■ \ j. 1 1 

2008 Oct 3 1 








V 


300 


1.5 


1.01-1.03 


— w \ / \ l ) /w j 1 

2008 Oct 3 1 








I 


600 


1.5 


1 no 1 f\ a 

1.02-1.04 


—\t\t\l) XT / \ O 

2008 Nov 02 


Scl2 S 


ni no c — i nn 

01:02:52.08 


00 c r\ i~\ a r\r\ 

-33:50:24.00 


B 


300 


1.3 


1 no 1 f\ a 

1.02-1.04 


"WWW) XT / \ O 

2008 Nov 02 








V 


300 


1.3 


1 no 1 no 

1.02-1.03 


~W W \ ( ) XT / \ O 

2008 Nov 02 








I 


600 


1.3 


1.01-1.03 


2008 Oct 3 1 


Scl3 N 


00:57:15.48 


00 01 ~\ t\^\ 

-33:32:22.92 


B 


300 


1.5 


1.13-1.18 


2008 Oct 3 1 








V 


300 


1.3 


1.15-1.20 


2008 Oct 3 1 








I 


600 


1.3 


1.16-1.22 


rs nAQ XT AO 

2008 Nov 02 


Scl3 S 


f\r\ rn \ c AO 

00:57:15.48 


01 ca in m 

-33:50:22.92 


B 


300 


1.8 


1.01-1.03 


~W W \ ( ) XT / \ O 

2008 Nov 02 








V 


300 


2.0 


1 no 1 f\ a 

1.02-1.04 


—\t\t\li XT nO 

2008 Nov 02 








I 


600 


1.7 


1 no 1 f\ a 

1.02-1.04 


2008 Oct 3 1 


Scl4 N 


n 1 n n /"v-~> n s~ 

01:00:03.96 


00 ci on /~ r\ 

-32:57:30.60 


B 


300 


1.5 


1.25-1.35 


— w \ / \ l ) /" \ .01 

2008 Oct 3 1 








V 


300 


1.3 


1.25-1.40 


2008 Oct 3 1 








I 


600 


1.5 


1.30-1.42 


-1 AAO XT n O 

2008 Nov 02 


Scl4 S 


n 1 nn no c\/~ 

01:00:03.96 


00 1 c on /~ r\ 

-33:15:30.60 


B 


300 


1.7 


1.05-1.08 


-\ AAO XT AO 

2008 Nov 02 








V 


300 


1.9 


1.06-1.09 


^» n xt no 

2008 Nov 02 








I 


600 


1.5 


1.07-1.10 


rs nAQ XT AO 

2008 Nov 02 


Scl5 S 


01:00:03.96 


O/i c on o a 

-34:25:30.36 


B 


300 


1.5 


1.17-1.23 


-1 AAO XT AO 

2008 Nov 02 








V 


300 


1.5 


1.18-1.27 


—\t\t\li XT / \ O 

2008 Nov 02 








I 


600 


1.4 


1.21-1.29 


rs nAQ XT AO 

2008 Nov 02 


Scl6 S 


01:05:40.20 


-33:50:03.12 


B 


300 


2.0 


1 no 10 

1.92-2.18 


~1 / W \ I ) XT / \ O 

2008 Nov 02 








V 


300 


1.8 


1 t \~i on 

1.97-2.30 


r\ AAO XT AO 

2008 Nov 02 








I 


600 


1.9 


m on 

2.07-2.39 


nAAO XT A 1 

2008 Nov 01 


Scl8 N 


ni Art no c\/~ 

01:00:03.96 


00 1 no 

-33:22:31.08 


B 


300 


1.5 


1.05-1.08 


rs nAQ XT m 

2008 Nov 01 








V 


300 


1.4 


1.04-1.07 


-\ / w w > XT / \ 1 

2008 Nov 01 








1 


600 


1.4 


1.03-1.06 


2008 Nov 01 


Scl9N 


00:57:16.92 


-32:57:23.40 


B 


300 


1.8 


1.00-1.01 


2008 Nov 01 








V 


300 


1.7 


1.00-1.01 


2008 Nov 01 








I 


600 


1.4 


1.00-1.01 


2008 Nov 01 


ScllON 


00:57:14.40 


-34:07:22.80 


B 


300 


1.6 


1.01-1.03 


2008 Nov 01 








V 


300 


1.5 


1.02-1.04 


2008 Nov 01 








I 


600 


1.4 


1.02-1.04 


2008 Nov 01 


Sclll N 


01:02:53.16 


-34:07:23.52 


B 


300 


1.5 


1.05-1.11 


2008 Nov 01 








V 


300 


1.5 


1.07-1.12 


2008 Nov 01 








I 


600 


1.3 


1.10-1.14 


2008 Nov 01 


Scll2N 


01:02:51.00 


-32:57:24.12 


B 


300 


1.3 


1.26-1.33 


2008 Nov 01 








V 


300 


1.4 


1.27-1.35 


2008 Nov 01 








I 


600 


1.4 


1.30-1.38 



Subsequently, we remove these bad pixels from the 
image by interpolating over them, preventing them 
from propagating into adjacent pixels during the re- 
sampling process. 

Due to the radial distortion in the large field of view 



the pixel scale is not constant across the field-of- 
view. Therefore the images had to be resampled to 
a common uniform pixel scale before stacking dif- 
ferent exposures that have been offset with respect 
to each other because of telescope pointing errors 
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Table 3: Photometric solution parameters and uncertainties. 



Band & colour 


4 


a 


P 


4m 


B,B-V 


5.133±0.009 


0.246±0.001 


-0.077+0.001 


I,V-I 


5.501±0.001 


0.062+0.001 


0.031±0.001 


V,B-V 


4.917±0.001 


0.127±0.001 


0.049+0.001 


V,V-I 


4.924±0.002 


0.129±0.001 


0.030+0.001 


0.9m 


B,B-V 


7.969±0.013 


0.309±0.006 


0.069+0.003 


I,V-I 


8.749±0.011 


0.031 ±0.005 


-0.028+0.002 


V,B-V 


7.936±0.013 


0.084±0.007 


-0.039+0.007 


V,V-I 


7.928±0.012 


0.084±0.007 


-0.046+0.006 



between multiple observations of the same field. 
This requires a precise World Coordinate System 
(WCS), which provides the rotation, scale and off- 
sets for mapping pixel coordinates to sky coordi- 
nates. The WCS header assigned to each image by 
the CTIO 4m-Blanco telescope was not sufficient 
to accurately align the images. Therefore, we ob- 
tained a new WCS solution using stars from the 
2MASS point source catalog dSkrutskie et al.ll2006l) 
that were found in ou r field-of-view. The solutio n 
uses a ZPN projection (Calabr etta & Greiserj2 002), 
which was fitted on each chip individually, while 
taking into account it's position with respect to a 
reference point (the center of the mosaic). 
From a comparison of the final calibrated positions 
with 2MASS positions we find that the WCS is ac- 
curate to +0.2" across the entire field of view. This 
level of accuracy (less than one pixel) is sufficient to 
align individual images. The images are resampled 
to a uniform pixel scale of 0.27" and then scaled to 
a common exposure time and sky level and stacked 
together to produce as deep a final image of each 
pointing as is possible. 

2.2. CTIO 0.9m 

The 0.9m telescope is equiped with a single 
2048x2046 pixel CCD detector, which has a field of 
view of 13.5' x 13.5' and a pixel scale of 0.396". We 
use standard bias and flat-field removal techniques 
in IRAF for the reduction. A predefined bad pixel 
mask provided by NOAO is used to flag bad pix- 
els, since the CCD's of the 0.9m telescope have sta- 
ble defects. A precise WCS is obtained, in order to 
match all stars observed in the 0.9m images to the 
same stars observed in the 4m images. As for the 4m 
data the WCS solution is derived using stars from 
the 2MASS catalog, using a ZPN projection of sky 



coordinates to pixel coordinates. The final solution 
is also accurate to +0.2" across the image, which 
corresponds to roughly half a pixel. This is suffi- 
ciently accurate to cross-correlate the 0.9m source 
catalogues with those of the 4m. 

3. Photometry 

Photometry is performed on all the re- 
duced i mages (0.9m and 4m) using 
DoPHOT dSchechter etail 1 19931) . using a ver- 
sion which is optimized to carry out accurate 
photometry down to the faintest levels. The 
method ology used in thi s section is described in 
detail in Saha et a □ (I2010h . 

DoPHOT uses an analytical function as a model 
PSF for describing different object types. Each 
model is characterized by six different parame- 
ters (X position, Y position, total intensity and 
three shape parameters), used both for the detection 
of objects and for determining their brightness. 
The initial shape parameters for star-like objects 
are determined separately for each image using a 
sample of bright, isolated stars. These parameters 
are updated after each consecutive pass through the 
image. 

The model PSF is not allowed to vary with position 
on the image. Instead, a fixed PSF is used, which is 
only allowed to vary from chip to chip. However, 
as the true PSF is expected to vary across the 
image (and also across each chip), we incorporate 
a position-dependency into the determination of the 
aperture corrections. Aperture magnitudes (which 
are independent of the PSF shape) are determined 
for bright isolated stars across each chip. The dif- 
ference between the PSF and aperture magnitudes 
is mapped as a function of position over the field 
of view, and this provides the correction factors for 
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Fig. 2: Top row: Calibrated PSF magnitudes plotted against photometric errors for the central pointing for B, V & 
I filters. Bottom row: Magnitude difference between stars in overlapping regions of different pointings, for B,V & I 
filters. 



DoPHOT PSF magnitudes. 

The fitting function used to map the difference 
between aperture and PSF magnitudes is different 
for 4m and 0.9m images. In the case of the 0.9m 
images, a simple quadratic surface is used, to allow 
variations across the image. In the case of the 
MOSAIC images a more complex surface is used, 
to allow variations across each chip to be treated 
independently, in addition to a quadratic surface 
across the whole image to fi t tilts or offs e ts be tween 
the individual chips (See ISaha et all d2010h for 
details). 

The photometry of the short and long exposures 
is combined by first cross-correlating the stellar 
positions. Subsequently, for stars that have are 
found in both exposures the star with the smallest 
photometric error is used in the combined catalog. 
Stars with bright magnitudes that were found only 
in the short exposure are also added to the catalog. 

3.1. Photometric Calibration 

We u se stand a rd sta r obse rvations of fields taken 
from lLandoltl (1 19921 l2007h with included Stetson 



standard stars (Stetson 2000) to obtain photomet- 
ric solutions for each night. The photometry is cal- 
ibrated to the Johnson-Cousins-Glass photometric 
system using a photometric solution of the form: 

magobs = mag true + A mag + a*x+P* colour 

where mag b s is the observed magnitude, mag true 
the true magnitude, x the observed airmass, colour 
the observed colour index used, A mag the zeropoint 
correction term, a the airmass coefficient and /3 the 
colour term. 

We determined the photometric solutions automat- 
ically with a Bayesian approach, which provides a 
statistically robust method to separate reliable mea- 
surements from inaccurate ones. We used a Monte 
Carlo method to determine the most probable pho- 
tometric solution, given the standard star measure- 
ments and taking into account the possibility of an 
unknown fraction of inaccurate measurements. 
Using this method, we obtain likelihood distribu- 
tions for the coefficients of the photometric calibra- 
tion, the fraction of outliers and the parameters of 
the outlier distribution. Furthermore, the likelihood 
of being an outlier is given for each standard star 
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measurement. The mean of the likelihood distribu- 
tion is adopted as the photometric solution, with the 
uncertainty of each coefficient given by the standard 
deviation of the corresponding distribution. The re- 
sulting photometric coefficients, along with their 
uncertainties are given in Table. [3] A detailed de- 
scription of the checks done to ensure a reliable pho- 
tometric calibration is given in Appendix 151 
The resulting catalogs of different bands are 
matched together to create one full catalog contain- 
ing 147456 stars in the V band down to 25.5 mag, 
72577 in the B band going down to 25.7 mag and 
116414 in the I band going down to 24.3. The top 
panels of Fig. [2] show the DoPHOT internal errors 
of all these stars as a function of their magnitude 
for the central pointing only. The other pointings 
display similar profiles, with the exception of the 
outer B band observations, which have shorter inte- 
gration times and thus go to different depths. At the 
bright end the best error is ^0.002, while at the faint 
end the error is ^0.2 for the faintest MSTO stars 
at the magnitude levels quoted above. The catalog 
contains only the stars for which all PSF parame- 
ters could be determined (type 1) and those that are 
found as mergers of two stars (type 3). 
The photometry of the outer pointings is placed on 
the same system as the central pointing by using 
stars in the overlapping region. The multiple mea- 
surements are used to fit an offset and colour term 
which put the photometry on the same system as the 
central pointing. Comparison of the magnitudes of 
stars which have two or more measurements in the 
overlapping regions are shown in the bottom panels 
of Fig. [2] The figure shows that the photometry of 
all pointings is calibrated to the same scale as the 
central pointing within the uncertainty of the pho- 
tometric measurements. The photometric accuracy 
of the final catalog coming from random errors is 
~0.002 mag for the brightest stars, going down to 
siQ.2 mag for the faint stars. 

3.2. Artificial star tests 

Given that our photometry only just extends beyond 
the magnitude of the oldest MSTOs it is important 
to quantify the completeness down to the faintest 
magnitudes. Incompleteness is caused by photomet- 
ric errors, crowding effects and the detection limit 
of the images, and results in the loss of stars at 
faint magnitudes. In order to derive the complete- 
ness fraction at different magnitudes we have per- 
formed artificial star tests on the deep stacked im- 
ages. Tests are done by inserting artificial stars in 



the images (a number of the order of 5% of the to- 
tal) and re-doing the photometry in the same way. 
The resulting catalog is matched to the observed 
and artificial input catalogs to obtain the fraction 
of stars recovered (and their measurement errors) at 
different magnitudes and colours. The completeness 
fraction for different filters is plotted for the central 
Sculptor pointing in Fig. [3] as a function of magni- 
tude. For the other pointings the completeness pro- 
file looks similar in the V and I filters. Due to bad 
seeing conditions the B completeness profile in the 
outer pointings drops off faster than shown in Fig. [3] 
For magnitudes brighter than shown in the figure 
the catalogs are assumed to be complete. Fig. [3] 
shows that the 50% completeness level is located at 
B«;25.2, V^24.8 and 1^24.0, which means our lim- 
iting magnitudes are about 1 magnitude deeper than 
the oldest MSTOs. 




1 ' ' ' ' ' ' =« : 3= 1 

18 19 20 21 22 23 24 25 26 27 
mag 



Fig. 3: Completenes fractions of the central Sculptor field 
as a function of magnitude for the B,V and I filter. For 
stars brighter then shown here the completeness factor is 
assumed to be 100%. 



3.3. The Colour-Magnitude Diagrams 

Using the full calibrated photometry catalog we can 
now produce Colour-Magnitude Diagrams (CMDs) 
to study different evolutionary features of the 
Sculptor dSph. Fig.|4]shows an (I,V-I) CMD of the 
inner part of Sculptor, where we have identified the 
different evolutionary features that can be seen in all 
the CMDs. 

Figs.|5]and|6]show calibrated (I,V-I) and (V,B-V) 
CMDs for the Sculptor dSph for different an- 
nuli of elliptical radius (r e //). A line indicat- 
ing the 50% completeness levels of these stars 
is also included. The elliptical radius is de- 
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Fig. 4: Sculptor (I,V-I) colour magnitude diagram of the central region overlaid with regions identifying the different 
stellar populations. BHB: Blue Horizontal Branch, RHB: Red Horizontal Branch, RGB: Red Giant Branch, BSS: Blue 
Straggler Stars, MSTO: Main Sequence Turn-Off. In the RGB box numerous AGB stars can also be observed. A 
box is also included which provides an estimate of the contamination by foreground stars from the Milky Way (FG). 
Photometric errors are shown as the size of the errorbar at the corresponding magnitude level. 



fined as the major axis of the ellipse centred on 
Sculptor (RA=01:00: 09, DEC=-33 :42:30 with el- 
lipticity e=0.32, see iMated dl998h V We assume 
a tidal radi us of 76.5 arcmin and core r adius of 
5.8 arcmin (llrwin & Hatz idimitrioul [19951) . No de- 
reddening is applied to the CMDs in Figs. [5] and |6j 
Instead, models and isochrones used to analyse the 
CMDs are reddened using the same extinction coef- 
ficient as obtained for the Sculptor dSph. 
The CMDs show only the stars for which all PSF 
shape parameters are determined. Therefore, the 
CMDs do not go as deep as the limits determined 
from the completeness profile (Fig. |3j. If we con- 
sider the stars for which not all the shape parameters 
have been determined as 'lost', the 50% complete- 
ness levels for this subsample of stars are B*23.9, 
V*23.5 and 1*22.8. We stress that this does not 
mean that 50% of the stars at these levels are not 
observed, but only that these are not considered in 
the subsequent analysis. For the (I,V— I) CMDs full 
coverage is obtained for radii out to r e //=0.8 deg and 
for the (V,B-V) CMDs out to r eU =0.6 deg. 



The CMDs display several features which remain 
clearly visible out to r c //=l deg (such as the RGB 
and BHB). The RHB, however, disappears around 
r c //=0.6 deg where the number of foreground stars 
becomes comparable to the number of RHB stars. A 
population of potential Blue Straggler Stars (BSS) 
remains visible out to r e //=l deg. The old MSTO 
stars remain visible at all radii, but are strongly af- 
fected by Milky Way foreground contamination be- 
yond r e //=l deg. The RGB has a spread in colour 
which is indicative of a spread in metallicity and/or 
age. A composite RGB bump is also present, for 
which a zoom in of the CMD is presented in Fig. [7] 
This is consistent with the presence of an extended 
period of star formation. 

To estimate the effect of foreground stars in 
the CMDs, we hav e used the Besancon mod- 
els (Robin et al. I2003h to construct a model Galactic 
population in the direction of the Sculptor dSph. A 
comparison of the number of foreground stars in the 
CMD (FG, see Fig.|Ui with the predicted foreground 
stars of the same colours from the Besancon mod- 
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Fig. 5: (I,V— I) Colour-Magnitude Diagrams of the Sculptor dSph for increasing elliptical radius. Error bars show the 
average photometric error at that level. The 50% completeness level for this subsample is indicated by the (red) solid 
line, the blue line shows a reference isochrone ([Fe/H]=-1.80,[a/Fe]=0.2, age = 11 Gyr) which runs through the 
middle of the RGB in the centre. 




Fig. 6: (V,B-V) Colour-Magnitude Diagrams of the Sculptor dSph for increasing elliptical radius. Error bars showing 
the average photometric error and a 50% completeness line (in red) are also included. A reference isochrone indicating 
the middle of the RGB in the centre is also shown as a blue line. 



els agree to within 10%, validating the use of the 
Besangon models to accurately predict the number 
of foreground stars anywhere in the CMD. 
Besides foreground stars a cloud of faint red stars is 



seen close to the detection limit in the (I, V— I) CMD, 
red ward of V-/ a: 1.0, which increases in number as 
we move out of the centre of Sculptor. These ob- 
jects are consistent with the integrated colours of 
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Fig. 7: Zoom in of the RGB bump region in the (V,B-V) 
CMD of the Sculptor dSph. A composite RGB bump is 
visible at 0.78<B-V<0.88 and 19.6<V<20.1. 



background spir al and irregular galax ies at a red- 
shift z=0.5-0.8 dFukugitaetal.lll995l) that remain 
unresolved due to a lack of spatial resolution and 
sensitivity. 

In order to disentangle the positions of the dif- 
ferent contaminating populations present in our 
CMDs, we have constructed a Colour-Colour dia- 
gram (Fig. [8^). The prominent evolutionary phases 
and foreground and background contamination have 
been colour-coded. The position of the same popu- 
lations are shown in the (I, V-I) CMD (Fig.[8t>), and 
the (V, B-V) CMD (Fig.lt). 
It is clear from Fig. [8] that the Milky Way fore- 
ground contaminates several features in the Sculptor 
CMDs, including the MSTO. The foreground stars 
become distinct from the main Sculptor popula- 
tion in Colour-Colour space for the reddest colours, 
as can be seen from the distinct 'plume' of stars 
in Fig. [8ji at (V-I)«1.6-3.4 and (B-V>1.2-1.8. 
The unresolved background galaxies are most dis- 
tinct from the Sculptor MSTO population in (V-I), 
making them show up as a separate 'cloud' of red 
stars, while in (B-V) they overlap with the MSTO. 
Therefore, having three filters helps in the interpre- 
tation of the CMDs. CMDs using different combi- 
nations of filters are also differently sensitive to age 
and metallicity effects (de Boeret al., 201 1, in prep). 

3.4. Structural Parameters 

Having quantified the most prominent sources of 
foreground and background contamination in our 
CMDs we can now investigate the spatial distri- 
butions of the different CMD features using star 
counts in foreground/background corrected regions. 
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Fig. 8: a) Sculptor (B-V), (V-I) colour-colour dia- 
gram (black points) overlaid with colour-coded CMD re- 
gions: Horizontal Branch (light-blue), RGB/SGB (green), 
MSTO (dark blue), MW foreground (red) and unresolved 
galaxy background (yellow), b) (V,B-V) and (I, V-I) (c) 
colour magnitude diagrams with the same colour-coding. 
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Fig. 9: Contour maps of the MSTO and BSS populations overlaid on the individual stars, shown as dots. The outer 
(red) ellipse is the tidal radius of the galaxy. The inner (red) ellipse is the largest ellipse not influenced by contour 
deformation caused by the proximity of the edge of observational coverage. The outermost contour in each plot is 
defined as 1% of the foreground corrected maximum intensity level. Subsequent contour levels increase by a factor 2 
for each contour. North is up, and East is to the left. 



We map the spatial distributions of different popu- 
lations in the Sculptor dSph using 2D iso-density 
contour maps. The stars are convolved with a nor- 
malized Gaussian with a width of <x=0.05 deg to 
generate intensity distributions. These distributions 
are then used to construct iso-density contours. We 
obtain the intensity level of the foreground stars by 
applying the same procedure to the Besancon model 
predictions, spread homogeneously across the field. 
Since the centre of the resulting contour maps 
was found to be off-centre with respect to the 
published RA and DEC of the Sculptor centre 
we decided to re-derive the structural parame- 
ters of Sculptor using our current data set. Using 
2D spatial Hess diagrams with a bin size of 
0.02x0.02 deg we investigated the radial variation 
of the central position, eccentricity and position 
angle with the IRAF task ELLIPSE. The anal- 
ysis reveals that the central position of Sculptor 
in our data set is located at RA=01:00:06.36, 
DEC=-33:42:12.6, which is west (-0.66'±0.06') 
and north (0.2 88'+0.05') of the central coordi- 
nates listed by iMateol (Q998). The eccentricity is 
found to be e=0.26+0.01 and the position angle 
PA=85.26±0.91, compare d to the literature val- 
ues o f e=0.32 and PA=99 (llrwin & Hatzidimitrioul 
1995). These parameters for our data set are adopted 
for the rest of the analysis in this paper. 
With these parameters the contour maps of the 
bright stars (down to the HB) confirm previous re- 



sults, showing that the outer contours are more ellip- 
tical than the inner and that the spatial distribution 
of the RHB and BHB populations are different. We 
show for the first time the spatial distribution maps 
of the MSTO and BSS population in Fig. |9j The 
BSS population looks similar in extent to the RGB 
and BHB. It shows the same change in ellipticity of 
the contours with different radii, as observed in the 
brighter features. 

A small overdensity is visible around £=-0.5, 
/7=0.1 in the MSTO stellar distribution in Fig. [9] 
which is more prominent in the spatial Hess di- 
agram (Fig. [Tol l. The photometric overdensity has 
about 30+5.5 stars per bin of 4.5x4.5 arcmin, com- 
pared to the mean number of stars of 13.6+3.7 
per bin at the corresponding elliptical radius. The 
feature coincides spatially with the stars of a ve- 
locity sub structure found in ESO/FLAMES spec- 
troscopy (lBattaglial f2007). which are overlaid on 
the Hess diagram in Fig.[l0] The velocity substruc- 
ture is made up of stars found in a cold structure at 
0.2 < T e n < 0.6 deg, with heliocentric velocities be- 
tween 128 and 142 km/s and a velocity dispersion of 
2.4+0.7 km/s. The spectroscopic metallicities cover 
a range of -2.8<[Fe/H]< -1.6. The substructure 
can be split up into two parts, one at 0.2 < r f « < 
0.3 deg (crosses in Fig. ITOb with a narrow metallic- 
ity range of [Fe/H]=- 1.96+0.05 and one at 0.3 < r eU 
< 0.6 deg (circles in Fig. [Tot with a broad metallic - 
ity range of -2.8<[Fe/H]< -1.6. The photometric 
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overdensity matches the spatial position of the ve- 
locity substructure at 0.3 < r e « < 0.6 deg. In terms 
of the CMD, the stars in the substructure match the 
same distribution as the main body of Sculptor stars, 
indicating that they are likely drawn from the same 
population, with a broad metallicity range. 
In order to investigate asymmetries in the spa- 
tial distributions of the different CMD populations 
we performed starcounts within annuli centred on 
Sculptor. For different r e // the fraction of stars 
West (£<0) of the centre versus East >0) of the 
centre is measured within an annulus. Fig.fTTIshows 
the fraction Nwestf^East f° r the full CMD (green 
circles), a selection of the CMD containing only 
Sculptor (black crosses), BSS population (blue as- 
terisks), foreground stars (cyan filled boxes) and the 
unresolved background galaxies (red plusses). 

MSTO spatial Hess diagram 




Fig. 10: Spatial Hess diagram of the MSTO population 
centred on Sculptor. Overlaid on the diagram are stars 
found in a velocity substructure (crosses indicate the sub- 
structure at 0.2 < i e u < 0.3 deg and circles the one at 0.3 < 
r e u < 0.6 deg.). A small photometric overdensity is visible 
at £=-0.5, f7=0.1, which coincides with one of the velocity 
substructures. 



The Sculptor populations (overall Sculptor, BSS) 
are almost symmetric for small elliptical radii, be- 
coming more asymmetric further out. The fore- 
ground stars show a symmetric distribution at all 
radii, consistent with Galactic population models. 
The unresolved background galaxies (see Fig. [8} 
display an anti-symmetric distribution with regards 
to the Sculptor features, which might be caused 
by obscuration of background galaxies by Sculptor 
stars, or a position-dependent morphological mis- 



classification of galaxies and stars. The asymmetry 
difference of BSS stars with regard s to the over- 
all Scu lptor population discussed by Mapelli et al. 
(2009) is not reproduced here and is likely due to the 
comparison with the full CMD at large radii instead 
of with the overall Sculptor population. The figure 
shows that the overall Sculptor population displays 
a slight asymmetry, with =sl0% more stars in the 
West than in the East, in the same direction as ob- 
served earlier by Eskridge (1988b). 
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Fig. 11: Fraction of Sculptor stars in the West (£<0) 
vs East (£>0) for different annuli of elliptical radius. 
Different populations are indicated with different colors, 
full CMD (green circles), overall Sculptor (black crosses), 
BSS (blue asterisks), foreground (cyan filled boxes) and 
background galaxies (red plusses) 



4. Interpretation 

Using the calibrated photometric catalog we can 
look at the behaviour of any stellar population over a 
large fraction of the galaxy. We can study the radial 
distributions of different CMD features and search 
for variations with elliptical radii. By analysing the 
CMDs of the MSTO region and combining this with 
information coming from spectroscopic surveys we 
can quantify the ages and metallicities of the differ- 
ent populations making up the Sculptor dSph. 



4.1. The Horizontal Branch 
We know from 



previous stu d ies of the Scu lptor 
" 19991: iHurlev-Kelleretal] 



dSph dMa iewski etal.1 
1999; Tolst ov et al.ll2004 that two distinct HB pop- 
ulations are present, with different radial distribu- 
tions. This can also be seen in our data in Fig. Q~2] 
showing the foreground corrected number of stars 
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in each bin of elliptical radius in the RHB and BHB 
populations. Since the number of foreground stars 
is not predicted to vary significantly with position 
across Sculptor the foreground correction is made 
by taking the predicted number of Galactic stars per 
unit area (as given by the Besan^on models) within 
the same magnitude and colour range as the CMD 
feature. This number is then multiplying by the area 
of the ellipses of different r e y, and subtracted from 
the observed star counts. The figure shows that the 
RHB population dominates in the central region, af- 
ter which it falls off more rapidly than the BHB dis- 
tribution. The BHB population is more uniformly 
distributed across the galaxy, and dominates only in 
the outer region. 




Fig. 12: Radial histograms of RHB (red, dashed) and 
BHB (blue, solid) populations in the Sculptor dSph. 



Besides the spatial distribution, the HB of Sculptor 
has also been extensively studied to look for vari- 
able stars. More than 90 RR Lyrae stars have 
been found in Sculptor, with average [Fe/H]«; 
—1.7 or lo wer along with 3 an omalous metal-poor 
Cepheids dKaluznv et alJll995l) . This confirms that 
the bulk of the stellar population in Sculptor is 
metal-poor. In addition, 2 more metal-rich short 
period variables (with [Fe/H]< -0.7) were found 
within the core radius, in accordance with the spa- 
tial distribution of the two HB populations. 

4.2. The Red Giant Branch 

The RGB of Sculptor is long known to have 
a larger colour spread than seen in Globular 
Clusters dDa Costall 19841) . consistent with it having 
experienced an extended period of star formation. 
However, no radial differences as clear as those on 
the HB have been found. This is because all popula- 
tions overlap on the RGB, especially at low metal- 



licity, and are thus not clearly distinguishable from 
each other. 

In Fig. [13] we show the RGB region in the 
(I,V-I) CMD for different elliptical radii. Overlaid 
on the CMDs are two isochrone s from the 
Dartm outh Stellar Evolution Database (Pott er et al.l 
2008). The isochrones are representative of the 
most metal-poor ([Fe/H]=-2.45) and most metal- 
rich ([Fe/H]=-0.90) stars making up the bulk of the 
population of Sculptor, as derived from RGB spec- 
troscopy dTolstov et al.ll2004t iBatta glial 120071) . The 
most metal poor stars found in Sculptor with these 
surveys have a metallicity of [Fe/H]< -2.5, which 
is not available in the Dartmouth Stellar Evolution 
Database. Instead, the most metal poor isochrone 
available in the library ([Fe/H]=-2.45) is used. The 
ages of the isochrones have been selected to fit the 
colours of RGB stars of corresponding metallicity. 
For increasing elliptical radius it can be seen that the 
brightest and reddest part of the RGB diminishes, as 
does the blue side of the Subgiant Branch. From the 
overlaid isochrones it can be seen that this is broadly 
consistent with the disappearance of the metal-rich 
component, as expected from the HB distributions. 
On the RGB, low resolution Ca n triplet (R-6500) 
spectroscopy is availa ble for &630 individual stars 
using VLT/FLAMES (Hbistov et alJ l2004llBattagUal 



120071: IStarkenburgetal.1120101) . We correlated this 
information with our photometric catalogue to com- 
bine the colors, magnitudes and [Fe/H] for a sam- 
ple of these stars. In Fig. [14] we show an (I,V— I) 
CMD of the RGB region for different radii over- 
laid with [Fe/H] metallicities for probable members 
stars to the Sculptor dSph where Ca n spectroscopy 
is available ((Starkenburg et al. 2010). Besides mem- 
ber RGB stars, numerous AGB stars are visible as 
well, along with a number of possible foreground 
stars. 

In the central region low resolution targets were se- 
lected only from stars that also had high resolution 
abundances, which explains the brighter magnitude 
limit in the top left panel of Fig. [14] Fig. [T4l shows 
that in the inner region stars are present with metal- 
licities ranging from [Fe/H]=-1.0 to [Fe/H]> -3.0. 
As we go outward in the galaxy the more metal rich 
components which are present on the bright, red part 
of the RGB in the inner 0.2 deg disappear. This is in 
agreement with the radial trend observed from the 
RGB photometry in Fig. [13] 
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Fig. 13: The distribution of stars on the Red Giant Branch for different elliptical radii. A metal-rich isochrone ([Fe/H]a; 
-1.10, [afFejx -0. 10, age=7Gyr, red) and metal-poor isochrone ([Fe/H]a -2.45, [a/Fe]s;0.40, age= 14Gyr, blue) from 
the Dartmouth Stellar Evolution Database have been overlaid in each panel. 



4.3. The Main Sequence Turn-Offs 

The MSTO region of the Sculptor dSph has previ- 
ously only been studied using small fields of view. 
Thes e studies determined an age ra nge of 13+2 
Gyr dDa Costal Il984l: iDolphinl |2002|) . We present 
here the first wide-field study of the MSTO popu- 
lation using homogeneous photometry that covers a 
large fraction of the entire galaxy. 
A close-up of the MSTO region in the (I, V-I) CMD 
for different elliptical radii is shown in Fig. Q3] 
overlaid with the same two isochrones used in 
Fig- H3 (one relatively metal poor and old and one 
metal rich and younger). For increasing distance 
from the centre we see that the MSTO proper- 
ties change. The inner regions display both young 
metal-rich stars and old metal-poor stars, while the 
outer regions lacks the young, metal-rich popula- 
tion. 

Using isochrones which match the spectroscopic 
metallicities on the RGB, we have selected regions 
of the MSTO corresponding to different metallici- 



ties and ages (see Fig. [ToT ). Regions were defined 
by selecting stars within a certain distance from 
each isochrone with the given age and metallic- 
ity. The range of [Fe/H] and [a/Fe] abundances of 
the isochrones were chosen to cover t he abundance 
trend derived from HR spectroscopy dTolstov et al.l 
120091 Hill et al., in prep), with an equal spac- 
ing from metal-rich ([Fe/H] =-1.00) to metal- 
poor ([Fe/H]=-2.20). The ages of the isochrones 
were determined by comparing them in the CMD 
to stars with the same abundance range. Fig. [16] 
shows the selected regions in the CMD as well as 
the resulting radial histograms along with the distri- 
butions of the RHB and BHB stars. From the MSTO 
radial distributions we see that the selected regions 
form a gradient from the metal-poor (BHB) to the 
metal-rich (RHB) distributions. This shows for the 
first time that this gradient is not only in metallic - 
ity, but also in age. The presence of an extended star 
formation history was expected from studies of the 
brighter CMD features, but now for the first time 
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Fig. 14: Sculptor CMD in the RGB region at different elliptical radii, r e u with Sculptor photometry as (black) dots, 
overlaid with larger (coloured) filled circles are stars for which spectroscopic me tallicities are available. The colours 
indicate the [Fe/H] metallicities of individual RGB stars (Starkenburg et al. 2010), and the scale is given at the side of 
each CMD. 



the effects of the age gradient can be seen in the ra- 
dial distribution. The oldest, most metal poor pop- 
ulations (with [Fe/H]«; -2.5) were formed roughly 
14 Gyr ago, while the youngest most metal rich pop- 
ulations (with [Fe/H]w -1.0) were formed approxi- 
mately 7 Gyr ago. 



4.4. Blue Straggler stars 

In the Sculptor CMDs (Fig.[5]&[6]> a BSS population 
appears to lie above the MSTO region (see labels 
in Fig. |4j». In terms of understanding the evolution- 
ary history of Sculptor it is important to understand 
what process produced these stars. Specifically, are 
they the result of star formation that continued un- 
til the much more recent times, after the main bulk 
of star formation 7-14 Gyr ago? There is a mech- 
anism of creating BSS stars proposed for globu- 
lar clusters res ulting from ma ss transfer between 



low m ass stars (McCrea 



cesses (Sigurdsson et al 



T%1 . or collisional pro- 
1994), which relates to an 



evolution of the old population in Sculptor. 
The radial distribution of the BSS population is 
shown in Fig. [T7J along with the distributions of 
RGB and MSTO. The distributions of the other 
features have been scaled to same total number of 
stars as the BSS distribution. There is a good agree- 
ment between the distributions of BSS and RGB or 
MSTO, which represent the overall population in 
Sculptor. This seems to indicate that the BSS popu- 
lation is not linked to any distinct age or metallicity, 
but likely also has the same age and metallicity gra- 
dient as seen in the other CMD features. 
No obvious trend of a radial gradient in the BSS 
population is evident in the CMDs with increas- 
ing radius, indicating that the BSS population can 
not easily be separated in different components. In 
Fig. [17] there is also no evidence of a central con- 
centration of the BSS population with respect to 
the overall Sculptor population, and hence no in- 
dications of collisio nal BSS formation as seen in 
Globular Clusters (iFerraro et al.l 119971) . This is in 
agreement with the fact that the central luminos- 
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Fig. 15: The distribution of stars on the Main Sequence Turn-Off for different elliptical radii. Two isochrones represen- 
tative of a metal-rich ([Fe/H]« -1.10, [a/Fe]a; -0.10, age=7Gyr, red) and a metal-poor ([Fe/H]« -2.45, [a/Fe]a;0.40, 
age=14Gyr, blue) population are overlaid in each panel. 



ity density of dS ph's (0.055 L /pc 3 for Sculptor, 
see Mateo (1998)) is significantly lower than ob- 
ser ved in globula r clusters (3200 L Q /pc~ 3 for M3, 
see lHarrisI (Il996l) ). This means the collision rates in 
dSph's must be much lower than in globular clus- 
ters, which makes the collisional BSS formation 
channel less likely to be dominant than the binary 
channel. 

In order to test if the BSS population is made up of a 
"normal" young population we modelled the distri- 
bution of BSS stars in the CMDs using isochrones, 
covering the metallicity range observed from spec- 
troscopy, assuming an age of 3 Gyr. Taking com- 
pleteness effects into account the observed popula- 
tion of «1600 BSS stars should lead to «35 stars on 
the upper RGB. Thus, from the photometry alone 
it is not possible to rule out that the BSS are a 
"normal" young population. However, if the BSS 
stars were a young population they would be ex- 
pected to be even more centrally concentrated than 
the older populations, if the trend of increasing cen- 



tral concentration with younger ages holds for these 
stars. Since this is not observed in Fig. [17] BSS 
formation through mass transfer seems most likely, 
which results in a distribution similar to that of other 
po pulations. This is consisten t with results found 
by iMapelli et all d2007l 120091) for several nearby 
dwarf spheroidal galaxies, including Sculptor. 



5. Discussion 

We have shown that the two distinct stellar popu- 
lations found in the Sculptor dSph by iTolstov et al.l 
(2004) using imaging of HB stars and spectroscopy 
of RGB stars are linked not only to a metallic- 
ity gradient but also to an age gradient going from 
roughly 7 to 14 Gyr ago. The Main Sequence Turn- 
Offs (Figs. [15] & [ToT ) show unequivocally that the 
outer regions of Sculptor lack significant numbers 
of young stars (<10 Gyr). This is supported by 
the spectroscopic Can triplet abundances (Fig. fT4l. 
which show a decrease in metal-rich stars with in- 
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creasing r e g. 

By using a combination of photometry and spec- 
troscopic abundance distributions, we determine the 
age gradient using the MSTO (Fig.fToTi. This shows 
that a radial gradient is present in age and metallic- 
ity, which is observed also as a bimodality in the HB 
distribution and as a spread on the RGB and MSTO. 
The importance of the MSTO is shown by the fact 
that only using this part of the CMD it is possible 
to constrain the age gradient accurately. Linking the 
age gradient to the existing metallicity gradient sug- 
gests that the observed change in HB morphology is 
du e not only to a metall icity gradient (as suggested 
bv lHarbeck et al.l (1200 ll) ). but that age is also an im- 
portant parameter. 

No obvious radial trend with age or metallicity is 
visible in the BSS population, but there are rela- 
tively few stars and there is no obvious way to sepa- 
rate them into distinct components. However, since 
its radial distribution does match the distribution 
seen in the overall Sculptor populations (Fig. [TTT i. 
it is possible that the gradient is also present there. 




0.2 0.4 0.6 0.8 

'ell ( de S) 



Fig. 17: Radial distributions of BSS stars (black crosses) 
compred to RGB (light blue filled circles), and 
MSTO (green filled squares). 

These results are consistent with the picture where 
the Sculptor dSph first formed an extended old 
metal-poor component, after which later genera- 
tions of progressively more metal-rich stars were 
formed ever more concentrated towards the centre 
of the galaxy. It seems that about 7 Gyr ago the gas 
supply ran out and the galaxy continued passively 
evolving up to the present time. 
The fact that the population gradient is still visi- 
ble today suggests that it is unlikely that the in- 
ner parts of Sculptor have suffered any significant 
tidal disruption, which could mix the different stel- 
lar populations together. Indeed, the spatial distri- 



butions (Fig. |9]l do not show any signs of obvious 
recent tidal interactions. This is in agreement with 
predic tions of theoretic al models of the Sculptor 
dSph (Sale s et aTll2010h . One aspect which remains 
unclear is w hat has caused t he increasing ellipticity 
with radius (lBattagliall2007l) . This may be due to a 
very early disruption, which would be visible only 
in the oldest outer parts of Sculptor, beyond radii 
studied here. Additionally, it could also be due to 
an encounter that only affected the outer regions of 
Sculptor. 

In a forthcoming paper the CMDs we have pre- 
sented here will be used to derive a detailed Star 
Formation History of the Sculptor dSph. This will 
be combined with high resolution spectroscopic in- 
formation (Hill et al., in prep) to derive the enrich- 
ment timescales of various elements, and how this 
fits into the star formation history timescale. 
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Fig. 16: On the left are selected regions on the MSTO in green for the labelled age and metallicity. The corresponding 
radial histograms are shown on the right (in green). Radial histograms of BHB (blue) and RHB (red) stars are also 
plotted. The radial distributions of the selected regions change from RHB-like for high metallicity to BHB-like for low 
metallicity. 
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Appendix A: List of standard field observations 

Observations of standard star fields were earned out during the observing runs in September 2008 and 
November 2009 in order to provide an accurate calibration o f the science observations. Fields were se- 
lected from the Landolt standard fields (lLandoltll2007l 1 19921) . to be used together with Stetson standard 
star photometry (IStetsonll2000h . Using the CTIO 4m/Blanco telescope observations of standard star fields 
were taken during the night in between science observations, in order to monitor the photometric fidelity 
throughout the night. Standard fields were selected to cover a broad range of airmass and colour, to allow 
adequate determination of the airmass coefficient and colour term in the photometric solution. In order 
to ensure of the most accurate calibration possible additional observations of standard stars and Sculptor 
fields were taken using the CTIO 0.9m telescope during 3 nights in 2008 under photometric conditions. 
Using these observations the calibration of the 4m Sculptor fields using the 4m standards can be checked 
and adjusted to the photometric solution obtained under photometric conditions on the 0.9m telescope. 
Table I A. ll lists the observations of standard fields taken with the 4m Blanco telescope, during the runs in 
2008 and 2009. The names of the standard fields are given, along with information on exposure times, 
seeing conditions (on image) and airmass. Table IAT21 lists the observations of standard fields taken with 
the 0.9m telescope, under photometric conditions during the 2008 run. Field names are listed, as well as 
exposure times, airmass and seeing conditions on image. 
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Table A.l: List of observed standard star fields with the 4m CTIO Blanco telescope during two observing runs. 
Information is given about exposure time, airmass and seeing conditions as determined on image. 



Field 


RA 


DEC 


Filter 


exp time 


seeing 


airmass 




J2000 


J2000 




sec 


it 




2008 


TPhe 


00:30:34.4 


-46:28:08 


B 


10, 8 


0.9-1.1,1.9 


1.05,1.08,1.20,2.66 








V 


5,4 


0.8-1.0,1.9 


1.05,1.08,1.17,2.62 








I 


4, 3 


0.8-1.0,1.8 


1.05,1.08,1.17,2.55 


MCT 0401-4017 


04:03:01.2 


-40:11:28 


B 


10,8 


0.9,1.3 


1.02, 1.03 








V 


5,4 


1.0-1.1 


1.02, 1.03 








I 


4,3 


1.0-1.1 


1.02, 1.03 


MCT 2019-4339 


20:22:51.5 


-43:29:24 


B 


10 


1.3 


1.1 








V 


5 


1.4 


1.1 








I 


4 


1.1 


1.1 


JL 82 


21:36:01.0 


-72:47:51 


B 


10 


1.6,1.8 


1.57, 1.91 








V 


5 


1.5,1.8 


1.56, 1.92 








I 


4 


1.2,1.5 


1.55, 1.93 


2009 


TPhe 


00:30:34.4 


-46:28:08 


B 


30, 10 


0.9-1.1 


1.04,1.05 








V 


30,5 


0.8-1.0 


1.04,1.05 








I 


30,4 


0.7-1.0 


1.04,1.05 


MCT 0401-4017 


04:03:01.2 


-40:11:28 


B 


10 


0.9,1.7 


1.42,1.44 








V 


5 


1.0,1.3 


1.41,1.44 








I 


4 


0.9,1.3 


1.40,1.45 


MCT 0550-4911 


05:51:53.9 


-49:10:31 


B 


10 


1.0,1.6 


1.09,1.73 








V 


5 


0.9,1.2 


1.10,1.72 








I 


4 


1.0,0.9 


1.10,1.70 


WD 0830-535 


08:31:56.5 


-53:40:52 


B 


10 


1.5 


1.13 








V 


5 


1.1 


1.13 








I 


4 


0.9 


1.12 


LB 1735 


04:31:34.7 


-53:37:13 


B 


10 


1.0 


1.1 








V 


5 


0.9 


1.1 








I 


4 


0.8 


1.1 


L92 


00:55:38.0 


+00:56:00.0 


B 


30,5 


1.5 


1.18 








V 


30,5 


1.2 


1.17 








I 


30,5 


1.3 


1.17 


L101 


09:56:39.0 


-00:19:53.0 


B 


30,5 


1.8 


1.42 








V 


30,5 


1.6 


1.45 








I 


30,5 


1.2 


1.49 


PG 2213-006 


22:16:33.5 


-00:23:34.7 


B 


30,5 


1.5 


1.27 








V 


30,5 


1.2 


1.29 








1 


30,5 


1.2 


1.30 


PG 0231+051 


02:33:55.3 


+05:14:22.9 


B 


30,5 


1.3 


1.60 








V 


30,5 


1.1 


1.57 








I 


30,5 


1.0 


1.54 


L95 


03:53:58.0 


+00:03:30.0 


B 


30,5 


1.5 


2.10 








V 


30,5 


1.4 


2.00 








I 


30,5 


1.4 


1.90 


L95 


06:51:58.7 


-00:24:18.7 


B 


30,5 


1.0 


1.17 








V 


30,5 


1.0 


1.17 








1 


30, 5 


0.9 


1.16 
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Table A. 2: List of observed standard star fields with the 0.9m CTIO telescope. Information is given about exposure 
time, airmass and seeing conditions as determined on image. 



rield 


KA 

J2000 


Uht 
J2000 


Filter 


exp time 

sec 


seeing 


airmass 


JL 163 


00:10:42.6 


-50:14:10 


B 


150 


1 A 1 O f\ 

1.4, 1.8-2.0 


1.16,2.31 








V 


150 


1.3, 1.8-2.1 


1.15, 2.39 








I 


300 


1.2, 1.6-1.7 


1.13, 2.53 


T Phe 


C\r\ o r\ o A A 

00:30:34.4 


-46:28:08 


B 


150 


1 A 1T -1 -1 A 

1.4-1.7, 2.2-2.4 


1 l\ A 1 t n 1^1 \ r>£l 1 n 

l. 04-1. 12, 1.21-1.26, 1.57 








V 


150 


1.3-1.7, 2.0 


1 (\A 1 1 O 1 in 1 OO 1 £LC\ 

1.04-1.13, 1.20-1.28, 1.60 








I 


300 


1.3-1.7, 2.2-2.4 


a f\ A i i r 1 1 n i a 1 s c 

1.04-L15, 1.18-L34, 1.65 


MCT 0401-4017 


04:03:01.2 


-40:11:28 


B 


150 


1.4-1.5, 2.1 


1.04-1.07, 1.16-1.24, 2.16 








V 


150 


1.3-1.6, 2.4 


1.05-1.08, L18-L20, 2.08 








I 


300 


1.2-1.5, 2.0 


1.04-1.09, 1.21-L18, 2.00 


LB 1735 


04:31:34.7 


-53:37:13 


B 


150 


1.3,2.3 


1.H-1.17, 1.38 








V 


150 


1.4, 1.9 


1.H-L18, 1.37 








I 


300 


1.2, 2.0 


1.12-L20, 1.34 


LSE 259 


16:53:45.1 


-56:00:16 


B 


150 


4.0 


2.16 








V 


150 


2.2 


2.21 








I 


300 


2.3 


2.29 


JL82 


21:36:01.0 


-72:47:51 


B 


150 


1.8 


1.38-1.40 








V 


150 


1.5-1.7 


1.39-1.41 








I 


300 


1.4 


1.40-1.42 


JL 117 


22:54:55.5 


-72:22:59 


B 


150 


1.5 


1.35 








V 


150 


1.4 


1.35 








I 


300 


2.0 


1.35 



23 



TJ.L. de Boer et al.: Deep Wide-Field Imaging of the Sculptor dSph 

Appendix B: Photometric calibration 

Obtaining an accurate photometric calibation is essential for a reliable interpretation of a CMD. For exam- 
ple the absolute luminosity of a Main Sequence Turn-off star correlates with age. The observing strategy 
was chosen in such a way as to combine deep 4m science images with 0.9m calibration images taken under 
photometric conditions. These images are used to ensure the most accurate photometric calibration possi- 
ble for the science images. 

A comparison of the calibrated magnitude and true magnitude of the standard stars used to obtain the pho- 
tometric solution versus observed magnitude (for both the 4m and 0.9m data) is shown in Fig. IB. II Similar 
plots showing the magnitude difference versus airmass and colour are shown in Figs. IB.2l and lB.3l respec- 
tively. A line indicating the zero level residual (black) and the mean of the residual (blue) is shown, along 
with errorbars denoting the average photometric error. These plots show the accuracy of the calibration 
applied to the science data. The standard deviation of the residuals is <x B «0.033, cry ^0.028, 07 ^0.032 
for the 4m data and cr B ^0.038, cry -0.028, 07 *0.029 for the 0.9m data. Fig. lB.ll shows that the zeropoint 
correction is consistent across all magnitudes used in the photometric solution, with the mean level of the 
residuals (blue line) being less than the average photometric error. Figs. IB.2l and lB.3l show that the effects 
of colour and airmass are well taken care of, with the mean level of the residuals (blue line) being less than 
the average photometric error. Residuals in the B band observations of the 0.9m telescope are most offset 
from zero, due to the poorer seeing conditions and larger photometric errors in the B band. However, the B 
band residuals are consistent with zero in all plots, given the average photometric error in B, showing that 
an accurate calibration was obtained. 

A second check of the photometry is done by comparing aperture and PSF magnitudes. We determined aper- 
ture magnitudes for small samples of stars in standard fields and in the central Sculptor pointing. A photo- 
metric solution was obtained from the aperture photometry, and used to calibrate the aperture magnitudes 
in the Sculptor central field. Then, a comparison of the PSF and aperture magnitudes is made (shown in 
Fig. lB.4l > for both the 0.9m and 4m data. The standard deviation of the residuals is cr B ~0.101, cry ~0.040, 
07 -0.033 for the 4m data and cr B «0. 108, cry ^0.028, 07 ^0.029 for the 0.9m data. Fig. IB11 also shows 
that the differences between aperture and PSF magnitudes are negligible within the average photometric 
errors, giving confidence to the reliability of our PSF photometry. The B band residual is once again most 
offset from zero (by 0.05), but consistent with zero given the photometric errors (+0.061). 
Next, the calibrated 4m and 0.9m PSF mags and aperture mags are compared to check whether the 4m 
data is consistently calibrated with the 0.9m data taken under photometric conditions. A comparison of 
both data sets for the central Sculptor field is shown in Fig. IB. 5 1 for PSF and aperture magnitudes re- 
spectively. Solid black errorbars have been overplotted for the PSF magnitudes, indicating the average 
error on the magnitude difference at bright, intermediate and faint magnitudes. For all three selections 
the mean of the residual is smaller than the error on the magnitude difference (from bright to faint: 
B: 0.0043<0.0617,0.0079<0.0880,0.0095<0.1546, V: 0.0099<0.0176,0.0154<0.0484,0.0187<0.0890, 1: 
0.0002<0.0125,0.0072<0.0272,0.0185<0.0844). For the aperture magnitudes the mean of the residuals (B: 
0.0170 V: 0.0321 1: 0.0236) is comparable to the average error (B: 0.0759 V: 0.0292 I: 0.0191). The rel- 
atively high residuals are possibly due to the low number of stars with aperture magnitudes used in the 
determination of the photometric solution for the 0.9m data. 

The above figures show that the 4m calibration properly takes into account zeropoint, airmass and colour 
effects. Furthermore, the calibration is consistent between aperture and PSF magnitudes, given the pho- 
tometric errors. The absolute calibration of the 4m data is consistent with that of the 0.9m data for PSF 
magnitudes and to a lesser extent for aperture magnitudes. Thus, an accurate absolute calibration is achieved 
for the PSF magnitudes used in the final photometry catalog. The final catalog has an average accuracy due 
to random errors for all filters of w +0.002 for the brightest stars, while at the faint end the accuracy is 
« ±0.2. The accuracy of the calibration varies for different magnitudes and filters, but the comparison of 
calibrated and true magnitudes shows an accuracy of ^0.04 mag or better across the magnitude range used. 
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Fig.B.l: Comparison between calibrated and true magnitudes for the standard stars used to obtain the photometric 
solution. Each plot shows magnitude residuals vs. true magnitude for different filters for the 4m data (left) and 0.9m 
data (right). Stars classified as outlyers when determining the photometric solution are marked in green. The average 
standard deviation of the residuals over all filters is aO.033 for the 4m data, and «0.03 for the 0.9m data. A line 
indicating the zero level residual (black) and the mean of the residual (blue) is also shown (visible only if there is an 
offset from zero), along with errorbars denoting the average photometric error. The zeropoint correction residuals are 
consistent with zero across all magnitudes, given the average photometric error. The mean of the residual in the 0.9m 
B (upper right) solution is most offset from zero (by 0.024) due to seeing conditions and large photometric error in 
B. However, the residuals are still within the average photometric error (±0.032), indicating an accurate calibration is 
achieved. 
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Fig.B.2: Magnitude residuals vs. airmass for photometric standard stars for the 4m data (left) and 0.9m data (right). 
See Fig. lB.ll for details. No trend with increasing airmass is visible, showing that the effects of airmass are well taken 
care of in the photometric calibration. 
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Fig. B.3: Magnitude residuals vs. colour for 4m (left) and 0.9m (right). standard stars. See Fig. lB.ll for details. No trend 
with colour is visible, showing that the photometric calibration correctly takes care of the effects of the colour term. 
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Fig. B.4: Comparison of aperture and PSF magnitudes in the central Sculptor pointing for the 4m data (left) and 0.9m 
data (right). See Fig. lB.ll for details. The residuals are consistent with zero, given the photometric errors, showing that 
the calibrations using aperture and PSF magnitudes agree. 
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Fig.B.5: Comparison between 4m and 0.9m photometry for the central Sculptor pointing for PSF (left) and aper- 
ture (right) magnitudes. A plume of short period RR Lyrae variable stars is visible in all bands in the PSF data at 
B~ 20.5 - 21; V~ 20 - 20.5; I~ 19.5 - 20.5. Solid black errorbars overplotted in the left-hand figures indicate the 
average error on the magnitude difference at the corresponding magnitude. The residuals of the PSF magnitudes are 
consistently smaller than the error on the magnitude difference, while for the aperture magnitudes the residuals are 
similar to the error. 
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